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ABSTRACT 

Elliptical galaxies and their groups having the largest L x /Lb lie close to the locus L x = 4.3 x 
10 43 (Lb/10 11 Lbq) 1 ' 75 expected for closed systems having baryon fractions equal to the cosmic mean 
value, fb ~ 0.16. The estimated baryon fractions for several of these galaxies/groups are also close 
to fb = 0.16 when the gas density is extrapolated to the virial radius. Evidently they are the least 
massive baryonically closed systems. Gas retention in these groups implies that non-gravitational 
heating cannot exceed about 1 keV per particle, consistent with the heating required to produce the 
deviation of groups from the L x — T correlation for more massive clusters. Isolated galaxies/groups 
with X-ray luminosities significantly lower than baryonically closed groups may have undermassive 
dark halos, overactive central AGNs, or higher star formation efficiencies. The virial mass and hot 
gas temperatures of nearly or completely closed groups correlate with the group X-ray luminosities 



and the optical luminosities of the group-centered elliptical galaxy, i.e. M v 



oc 



1.33 



an expected 



consequence of their merging history. The ratio of halo mass to the mass of the central galaxy for 
X-ray luminous galaxy/groups is M V i r /M* ~ 80. 

Subject headings: galaxies: elliptical and lenticular, CD - X-rays: galaxies - galaxies: clusters: general 
- X-rays: galaxies: clusters - galaxies: cooling flows 
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1. INTRODUCTION 

Massive elliptical galaxies with similar optical lumi- 
nosities have hot gas X-ray luminosities that range over 
two orders of magnitude. The origin of this scatter, 
shown in Figure 1 , has received much attention but a full 
understanding remains elusive. There is evidence that 
gas loss by ram pressure (and tidal) stripping has reduced 
L x /Lb in elliptical galaxies or groups orbiting within 
rich clusters of galaxies (Biller et al. 2004; Machacek 
et al. 2005; Sun et al. 2005). However, an enormous 
range in L x / Lb also prevails among non-interacting el- 
lipticals that are isolated or at the centers of isolated 
galaxy groups. The correlation between the spatial ex- 
tent of the X-ray emission and L x /Lb suggests that the 
driver for this scatter is a variation in the virial mass 
M V i r of the halo that surrounds otherwise similar ellip- 
tical galaxies (Mathews & Brighenti 1998). The virial 
mass M v i r and radius r„j r are found by fitting dark NFW 
halos to the total mass distribution derived from X-ray 
observations of the hot gas density and temperature in 
10 kpc < r < r V i r , assuming hydrostatic equilibrium. 

To gain further insight into the broad range of X-ray 
emission from optically similar galaxies, we draw at- 
tention here to those ellipticals with the largest X-ray 
luminosities. These isolated galaxy/groups have been 
variously referred to as "[X-ray] overluminous elliptical 
galaxies (OLEGS)" (Vikhlinin et al. 1999) or "fossil 
groups" (Ponman et al. 1994). The concept of fossilized 
groups is meant to imply that they are relics of merg- 
ing among galaxies in a group environment, although 
all elliptical galaxies may qualify for this designation. 
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Jones et al. (2003) provide an operational definition 
for fossil groups in terms of the magnitude difference 
between the first and second brightest group galaxies. 
For our purposes here we simply consider those elliptical 
galaxies with the largest L x /Lb in the L x — Lb plot, 
many of which have been previously regarded as fossils 
or OLEGS. We then note that several of the best studied 
of these galaxies have nearly the same baryon mass frac- 
tion as the most massive galaxy clusters and the WMAP 
value, fb = 0.16 (Spergel et al. 2003), i.e. they appear 
to be baryonically closed. Most baryons are in the hot 
intragroup gas. 

The data in Figure 1 are mostly taken from O 'Sullivan 
et al. (2001) (open squares), but we have added ad- 
ditional X-ray luminous ellipticals assembled from more 
recent observations (filled symbols) with properties listed 
in Table 1. These X-ray luminous systems define the up- 
per envelope of the luminosity distribution in the L x — Lb 
plane. While all estimates of the baryon mass fraction 
fb require uncertain extrapolations beyond the observa- 
tions to the virial radius r V i r , fb for several X-ray lumi- 
nous groups in Table 1 indicate near or complete baryon 
closure. All data have been scaled to Hq = 70 km s _1 
Mpc -1 . 

2. NON-GRAVITATIONAL HEATING IN THE NGC 5044 
GROUP 

Galaxy groups and poor clusters with M V i r <3x 10 14 
Mq and kT < 4 keV are known to deviate systematically 
below the L x — T relation established by more massive 
clusters, suggesting additional non-gravitational energy 
by cosmic preheating or AGN activity (as reviewed by 
Voit 2005). Consequently, it is remarkable that groups 
in Table 1 with kT ~ 1 — 2 keV have survived with most 
or all of their baryonic component intact. 

NGC 5044 is a good example of such a group. For 
the purpose of this discussion, we have made a prelim- 
inary mass model of NGC 5044 based on gas density 
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and temperature profiles observed to ~ 0.3r„j r « 300 
kpc from Buotc et al. (2003, 2004, & 2006 in prep.). 
In the central regions the azimuthally averaged gas den- 
sity n et0 bs was replaced with n e j = n e ^ obs / f 1 ^ 2 where 
f(r) = 1- 0.642 exp(-r fepc /20) (Buote et'al. 2003) is the 
filling factor of the denser gas component at each radius 
responsible for most of the observed emission. The model 
was constructed by first setting the stellar parameters - 
a de Vaucouleurs profile with luminosity Lb — 4.5 x I0 10 
Lbq, effective radiius R e = 10 kpc and stellar mass to 
light ratio Tb — 7.5 - that establish the total stellar mass 
M*e — 3.4 x I0 11 Mq and potential. The dark halo is as- 
sumed to have an NFW mass profile with an adjustable 
virial mass M vir and concentration c = 433M~ r 125 ex- 
pected for this mass (Bullock et al. 2001). The equation 
of hydrostatic equilibrium is integrated for n e (r) , fixing 
the gas temperature T(r) to fit observations and extrap- 
olating to larger radii in a logT — logr plot. M vir and 
the innermost gas density are varied until an excellent fit 
is achieved to the n e j(r) profile throughout the observed 
region. The resulting virial mass, M v i r = 4.0 x 10 13 Mq, 
is similar to our previous estimate (Buote et al. 2004) 
and the virial radius r vir = (3M t , il ./47rAp c ) 1 / 3 = 841 kpc 
with A = 104 and p c = 9.24 x 10~ 30 gm cm" 3 . When the 
observed gas density profile in NGC 5044 is extrapolated 
to r vir (Buote et al. 2004; 2006 in prep.), maintaining the 
same power law n e — 0.66r^ c 45 cm~ 3 observed in the re- 
gion 100 < r < 300 kpc, we find that the total gas mass is 
M g = 4.8 x 10 12 Mq, in agreement with the mass model. 
The mass fraction in gas is f g « 0.11. This corresponds 
to a baryon ratio f b w / s /(l — 0.12) = 0.13, assuming 
a (conservative) star formation efficiency of 12% (Lin & 
Mohr 2004). At least 80% of the initial baryons in NGC 
5044 is still bound to the group. Evidently, the non- 
gravitational heating received by the gas is < 20% of the 
gas binding energy, E bind — 9.6 x 10 61 ergs. 4 For sim- 
plicity we assume that the percentage difference between 
the observed f g and the value f g w (1 — 0.12)/b = 0.14 
expected from WMAP is proportional to the amount of 
non-gravitational energy that the gas received as a per- 
centage of E bind . 

The gas heating efficiency associated with accretion 
onto the central black hole in NGC 5044 must be con- 
sistent with gas retention. The central galaxy with mass 
M* = 3.4 x 10 11 Mq is expected to contain a black hole 
of mass M bh « 7.6 x lO^M 1 - 12 = 6.2 x 10 8 M Q (Haering 
& Rix 2004) . During the accretion history of the central 
black hole, suppose that a fraction rjh of the rest energy 
Mft^c 2 heats the intragroup gas, then gas retention at the 
80% level suggests that rj h < 0.2E bmd /M bh c 2 = 0.016, al- 
though some of this energy will be radiated away. The 
energy radiated by NGC 5044 during several Gyrs is 
E ra d = L x t = 4.3 x 10 60 (i/5 Gyrs) ergs and, since no gas 
is observed to cool below ~ T vir /3 in NGC 5044 (Buotc 
et al. 2003), the minimum accretion heating efficiency 
is Vrad ~ E rad /M bh c 2 ~ 0.004(t/5 Gyrs). Evidently 
only a tiny fraction of the accretion energy released by 
the central black hole (~ Q.lM^c 2 ) can have heated the 
intragroup gas in NGC 5044. Nevertheless, substantial 

4 The binding energy is found by computing the double inte- 
gral E bind = ^ vir p g Airr 2 dr J r °° g(r')dr' where p g and g(r) = 
GM(r)/r 2 are the gas density and gravitational acceleration from 
our mass model for NGC 5044. 



ongoing AGN-related heating is currently observed near 
the center of the NGC 5044 group (Buote et al. 2003; 
Mathews, Brighenti & Buote 2004). 

The non-gravitational energy received by the hot intr- 
acluster gas from supernovae can be estimated from the 
total mass of iron observed in the NGC 5044 group, 

M Fe = 0.71z F e,ga S M gas + 0.71z Fe ^aM* (1) 

where the estimated total gas mass is M gas = 4.8 x 10 12 
Mq, 0.71 is the ratio of hydrogen mass to total mass 
including helium, and the total stellar mass in the group 
is g times larger than M*. The mass- weighted gas iron 
abundance in the NGC 5044 group is z Fe ^ gas = 0.16zi?e© 
and we adopt a mean stellar abundance z Fe .* = 0.5z Fe Q 
and z FeQ = 1.83 x 10~ 3 (Grevesse & Sauval 1998). The 
total iron mass can be expressed in terms of supernova 
yields (yu ~ 0.1 Mq and yi a ~ 0.7 Mq) and r\, the 
number of supernovae per Mq of initial stars formed, 

M Fe = M n {rniy u + i] Ia y Ia ). (2) 

We assume that none of the supernova iron either cooled 
(Brighenti & Mathews 2005) or was buoyantly expelled 
to ~ r vir . If all stars formed at high redshift with 
a Salpeter IMF between 0.08 and 100 Mq the num- 
ber of stars with mass above 8M© that become Type 
II supernovae is r\u = 0.0068 per Mq. If crM» is the 
current stellar mass in the group, the initial mass is 
M*j w crM»/(l - (3) where w 0.3 is the fraction of 
the original stellar mass ejected from stars (Brighenti & 
Mathews 1999). If the iron mass M Fe is eliminated be- 
tween the two equations above, we find that the num- 
ber of Type la supernova per Mq of initial stars is 
rji a w 6.7 x 10 _4 (T _1 . The total energy released by both 
types of supernovae is 

E sn = M„(rj H + Vla )10 51 » 1.1 x 10 61 ergs (3) 

where we assume 10 51 ergs per supernova and a = 3. 
Since E sn / E b i nd w 0.11 it is possible that supernovae 
energy could eject ~ 10% of the baryons from the group. 
But our assumption that all of the SNII energy is com- 
municated to the hot gas may not be plausible since it 
makes no allowance for radiation losses in SNII remnants. 
Of course a flatter IMF (e.g. Brighenti & Mathews 1999; 
Nagashima et al. 2005) could generate enough SNII en- 
ergy to eject 20% of the gas (e.g. Brighenti & Mathews 
2001). Clearly, the large baryon fraction f b observed in 
NGC 5044 and other X-ray luminous groups imposes a 
significant constraint on non-gravitational heating. 

3. L B AND THE VIRIAL MASS 

The X-ray luminosity of the groups with filled symbols 
in Figure 1 along the upper envelope of the L x —Lb distri- 
bution correlates with Lb as L x w 1.7 x 10 21 (Lb/£_bo) 2 
erg s _1 (also see Jones et al. 2003). If these groups are 
essentially baryonically closed, as we propose here, then 
Lb for the group-centered E galaxy should also increase 
with the virial mass. This correlation is shown in Fig- 
ure 2a where we plot those groups from Table 1 having 
known estimated M vir . Similar Lb — M vir correlations 
have been found from the 2MASS survey (Lin & Mohr 
2004) and the galaxy-galaxy lensing data of Cooray & 
Milosavljevic (2005) who find M vir oc £^ 33 , which agrees 
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well with the group data in Figure 2a 5 . The Lb ~ M V i r 
correlation arises because group-centered elliptical galax- 
ies grow by mergers as massive satellite group galaxies 
undergo dynamical friction in the dark halos. Progres- 
sively more massive halos contain more mergeable galax- 
ies so the final Lb of the group-centered elliptical in- 
creases with M v i r until M v { r > 10 14 Mq where dynam- 
ical friction is reduced by the small galaxy/halo mass 
ratio. All baryonically closed groups have been dynami- 
cally processed in this way. 

The ratio of the total halo mass to that of the central 
galaxy Myir/M* is of particular interest. The correla- 
tion M mr = 1.12 x 1O 14 (Lb/1O 11 L B0 ) :l - 33 M q shown in 
Figure 2a (dashed line), when combined with the stellar 
mass to light ratio mJ/L b = 7.1( J L i3 /10 10 ) ' 29 (Trujillo 
et al. 2004), results in M vir /M* « 80(M»/10 12 M Q ) - 03 . 
This ratio is larger than the dynamical masses deter- 
mined from the Sloan Digital Sky Survey: M v i r /M* ~ 
15 - 20 (based on M V „/L B from Prada et al. 2003) 
and M mr /M* - 7 - 30 (Padmanabhan et al. 2004) for 
E galaxies comparable to those in Table 1. This sug- 
gests that optically similar elliptical galaxies with larger 
L x may have somewhat more massive dark halos. The 
mean gas temperature for these groups plotted in Figure 
2b also correlates with optical luminosity of the group- 
centered galaxy as indicated by the dotted line with slope 
T oc L° B 60 . Since L x cx L% from the Table 1 data 
in Figure 1, both M„j r and T also correlate with L x , 
M v „ oc L° x 7 and L x cx T 33 . 

The X-ray luminosity of baryonically closed groups 
(M gas oc M vir ) scales approximately 
(n e )M gas oc M 2 ir /r 3 ir cx M vir [For the mean gas tem- 
peratures of the groups in Table 1 , 1 < kT < 3 keV, the 
bolometric X-ray emissivity is insensitive to temperature 
(Sutherland & Dopita 1993)] . However, the L x — M vir 
relation also depends on the NFW concentration c = 
450(M vir /M Q )-°- 128 (Bullock et al. 2001). We estimate 
L x (c,M vir ) for closed groups (f g = f b - /* « 0.14) by 
filling NFW potentials with isothermal gas with T — 
2.22 x 10 7 (M OTr /10 14 M Q ) - 54 K taken from Shimizu et 
al. (2003). The bolometric X-ray luminosity within rw 
is L x = 4.9 x 10 43 (M t ,„/10 14 M Q ) 1 - 3 erg s" 1 for 10 135 < 
M V i r 10 14 5 M Q , where we assume M V i r oc L B 33 from 
Figure 2a (dashed line). This locus of maximum X-ray 
luminosity, L x = 4.3 x 10 43 ( J L B /10 11 is©) 1 ' 75 erg s"\ is 
shown with a dotted line in Figure 1. If the gas temper- 
ature T cx T vir oc M 2 / r 3 then L x oc M 1 ^ 3 suggests that 
L x oc T 2 similar to Kaiser (1986). 

4. CONCLUSIONS 

We propose that galaxy groups lying near the upper en- 
velope of the L x — Lb distribution in Figure 1 are nearly 
or completely baryonically closed boxes similar to more 
massive clusters. This conclusion is supported by the 
baryon fraction estimates listed in Table 1 and the prox- 
imity of the observations in Figure 1 to the approximate 
dotted line locus for the maximum L x expected from 
baryonically closed groups observed to r vir . The pro- 
jected X-ray luminosity of our mass model for NGC 5044 
beyond ~ 0.1r„j r varies as L x (rproj ) oc t^' 2 ^ so we expect 



that L x for galaxies with filled symbol points in Figure 1 
will creep upward toward the dotted line when observed 
with larger apertures and more sensitive detectors. Nev- 
ertheless, we do not expect to find galaxy/groups in the 
future that lie significantly above the filled circles and 
squares in Figure 1. 

Baryonically closed groups provide interesting con- 
straints on the amount of non-gravitational heating ac- 
quired by the intragroup gas. To retain most or all of 
the gas in these groups, the gas heating by the central 
black hole (AGN) must be r] h M bh c 2 < 0.016M bh c 2 or < 1 
keV per particle, consistent with typical values ~ 0.7 — 1 
keV per particle required to account for deviations from 
self-similarity in the L x — T plot for clusters (e.g. Tor- 
natore et al. 2003; Voit 2005). We also find that the 
combined energy of all past supernovae is insufficient to 
remove significant amounts of intragroup gas unless the 
IMF is flatter than Salpeter (Brighenti & Mathews 1999, 
2001). 

Another necessary attribute of baryonically closed 
groups is that they are spatially isolated, i.e. they have 
not lost mass by ram-stripping during mergers with com- 
parable or larger systems. Spatially isolated E galax- 
ies and groups are of particular interest because of the 
strong limits they impose on non-gravitational heating. 
It is therefore remarkable that some isolated E galax- 
ies have much lower L x than ellipticals in baryonically 
closed groups. For example, in Figure 1 we mark with + 
symbols two isolated E galaxies found recently by Reda 
et al. (2004) that are near the bottom of the distribu- 
tion. This large variation of L x for isolated Es of similar 
Lb in Figure 1 may result from normal cosmic variance. 
It would be interesting to determine if these and other 
isolated galaxies have unusually undermassive dark ha- 
los (allowing winds), if they contain more energetically 
active (and massive) black holes or if their star forma- 
tion efficiencies are unusually large, since such variations 
could help explain why these isolated galaxy/groups are 
not baryonically closed. By this means it will eventually 
be possible to determine if the non-gravitational heating 
arises primarily from the central black hole. It would also 
be worthwhile to assemble M V i ri (T) and optical lumi- 
nosities for all isolated E galaxies and groups throughout 
the L x — Lb plane. 

Finally, we have shown that baryonically closed groups 
can inform us about the important relationship between 
the optical luminosity and mass of the group-centered 
galaxy and the mass of the surrounding (group) dark 
halo. The preliminary data currently available suggest 
that the dark halos are about ~ 80 times more massive 
than the central (non-cD) elliptical galaxy. These con- 
clusions can be explored further in the L x — Lb plane by 
considering all elliptical galaxies for which X-ray obser- 
vations provide accurate M v i r . 

Studies of the evolution of hot gas in elliptical galax- 
ies at UC Santa Cruz are supported by NASA grants 
NAG 5-8409 & ATP02-0122-0079 and NSF grants AST- 
9802994 & AST-0098351 for which we are very grateful. 
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5 M vir for NGC 6482 (the leftmost point in Figure 2a) may be influence of the dark halo mass and the much smaller stellar mass, 

significantly underestimated since the X-ray observations extend 
only to ~ 30 kpc, close to the expected transition between the 
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TABLE 1 
LUMINOUS GALAXY GROUPS 3 - 



galaxy 




log L Xt bol 


ttl 


log M v i r 


Jb 


rtct. 


group 




(erg s _1 ) 


(keV) 


( M Q> 






NGC5044 


10.76 


42.80 


1.2 


13.60 


~ 0.13 


1 


RXJ1159 


11.09 


43.05 


2.2 


14.15 


~ 0.04 


2 


WJ943.7 


10.94 


43.07 


1.7 


13.74 


~ 0.15 c 


3 


RXJ0419 


10.64 


42.75 


1.4 


~13.6 


~ 0.03 c 


4 


NGC6482 


10.47 


42.04 


0.7 


12.60 


~ 0.18 c 


5 


ESO3060170 


11.30 


43.82 


2.7 


14.25 


NA 


6 


RXJ1340 


11.02 


43.11 


2.3 


14.33 d 


NA 


2,10 


RXJ2114 


10.96 


43.01 


2.1 


14.27 d 


NA 


2 


RXJ2247 


11.13 


43.32 


2.8 


14.45 d 


NA 


2 


NGC1132 


10.73 


42.71 


1.0 


13.52 


NA 


7,8 


RXJ1416 


11.24 


44.05 


1.5 


NA 


NA 


9 


RXJ1119 


10.36 


41.94 


NA 


NA 


NA 


9 


RXJ1256 


11.16 


43.49 


NA 


NA 


NA 


9 


RXJ1331 


10.68 


42.48 


NA 


NA 


NA 


9 


RXJ1552 


11.12 


43.51 


NA 


NA 


NA 


9 


RXJ0116 


10.76 


42.94 


NA 


NA 


NA 


9 



a All data scaled to H = 70 km s _1 Mpc -1 

b References: 1 Buotc ct al. (2004); 2 Vikhlinin, A. ct al. (1999); 3 Rasmussen, J. & Ponman, T. J. (2004); 4 Kawaharada, M. et al. 
(2003); 5 Khosroshahi, H. G., Jones, L. R. & Ponman, (2004); 6 Sun, M. et al. (2004); 7 Mulchacy, J. S. & Zabludoff (1999); 8 Gastaldello, 
F. ct al. (2005); 9 Jones, L. R. et al. (2003); 10 Jones, L. R. ct al. (2001) 

c Estimated from gas mass ratio by scaling up by 1.11 

d Rcfcrcncc 2 provided values of the total mass M(r x ) within radius r x . We used the definition of r v i r (M v i r ) with c(M„; r ) and y x = 
cr x /r v i r to determine M v i r = M(r x )f(c)/f(y x ) from the NFW mass profile where f(x) = ln(l + x) — x/(l + x). 
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L °g L B ( L B,©) 



Fig. 1. — Plot of the bolometric X-ray luminosity and B-band optical luminosity for elliptical galaxies (RC3 type T < —4) from 
O'Sullivan et al. (2001) (open squares). Only X-ray detected galaxies are shown. The filled circles (squares) are X-ray luminous ellipticals 
with estimated (unknown) virial masses as listed in Table 1. Two isolated elliptical galaxies (NGC 3557 and NGC 4697) arc marked with 
+ symbols. The dashed line approximately represents the stellar X-ray emission from binary stars (Kim & Fabbiano 2004). The dotted 
line L x = 4.3 X 10 43 (L s /10 n L BQ ) 175 is the locus of maximum L x for NFW halos maximally filled with gaseous baryons (f g = 0.14); 
the error bars show the effect of a 1<t change in concentration c(M v i r ) expected from cosmic variation (Bullock et al. 2001). On average 
the filled symbols lie ~ 0.33 below the dotted line; since L x oc /J and the observed L x < L x (r v i r ), they have gas filling factors f g > 0.08. 
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Fig. 2. — Plots of optical luminosity of baryonically closed groups against virial mass and mean gas temperature, (a) The dashed line 
in the upper panel shows the variation M„; r oc L^ 33 based on weak lensing. (b) The dotted line shows the correlation T tx L 1 ^ 60 



